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We propose an interaction flow scheme that sums up the perturbation expansion of many-particle 
systems by successively increasing the interaction strength. It combines the unbiasedness of renor- 
malization group methods with the simplicity of straight-forward perturbation theory. Applying 
the scheme to fermions in one dimension and to the two-dimensional Hubbard model we find that 
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^v^j . at one-loop level and low temperatures there is ample agreement with previous one-loop renormal- 

ization group approaches. We furthermore present results for the momentum-dependence of spin, 
^ ' charge and pairing interactions in the two-dimensional Hubbard model. 
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I. INTRODUCTION 



1} ' The competition between different Fermi surface instabilities is a general feature of low-dimensional interacting 
^ , lattice electron systems. In one spatial dimension a number of analytical and numerical methods are available, 
whereas in two and higher dimensions perturbative renormalization group methods are often considered the least biased 
! t | approach to study such systems at weak coupling and low energy scales. In recent years functional renormalization 
group methodsii2i2i£i£i§iL£i2ii£ have been devised, capable of capturing the interplay of various low-energy scattering 
processes within certain approximations. They provide a basis for a qualitative comparison of the strengths of potential 
instabilities. For those RG schemes which are derived from an exact flow equationSdiiiiiSii^, the strategy is to reduce 
the degree of approximation step by step in order to achieve quantitatively more accurate results^ However, the 
practical experience along these lines shows that each scheme has pros and cons. While easily applicable, the one- 
loop momentum shell schemes have difficulties describing ferromagnetic Stoner instabilities. 10 Schemes with frequency 
!— 1 ' cutoff obscure the analytical structure of correlation functions, and the temperature-flow schemed needs to be adapted 
to the specific type of the interaction. Thus, as long as an all-purpose weapon is not available, it is desirable to have 
a larger number of distinct approaches to tackle the problem from different angles. 
. Functional RG methods are extremely flexible concerning the choice of the flow parameter entering the quadratic 
' part of the fermionic action. The cutoff RG schemes are constructed by multiplying this quadratic part with a scale- 
\Q , dependent cutoff function. The derivative of this function with respect to the scale defines momentum or frequency 
C*~) 1 shells, which are integrated out in the course of the flow. The temperature-flow RG instead uses the temperature as 
the flow parameter, and the flow describes the change of the interactions as the temperature is lowered. 

In this work we use the interaction strength as the flow parameter. It is instructive to switch on the interaction 
strength adiabatically and to follow the change in the correlations as the interaction is increased. This defines an 
interaction flow (IF). Straight-forward perturbation theory could in principle do this to any order, but casting the 

S idea in the form of the functional RG has the compelling advantage that all one-loop corrections are summed up 
consistently. Therefore the IF scheme combines the unbiasedness of one-loop RG methods with a straight-forward 
■"^J ■ relation to perturbation theory. 

In this paper we apply the IF method to fermions in one dimension and to the two-dimensional (2D) Hubbard 
model. In the latter numerous instabilities may compete, and we can compare the nature of the infrared instabilities 
suggested by IF with results obtained from other RG approaches. Furthermore, the 2D Hubbard model is of general 
interest, in particular because of its potential relation to high-T c and other unconventional superconductors. However, 
we emphasize that the IF method is rather general and should be widely applicable to many kinds of interacting many- 
particle systems. 

The paper is organized as follows: In Sec. [H]we present the formalism for the IF scheme. In Sec. IHII we apply it to 
interacting fermions in one spatial dimension. In Sec. II VI we describe the implementation in two spatial dimensions 
for the Hubbard model on the square lattice and compare the results to those obtained in other one-loop RG schemes. 
Further more, we show new results for the detailed wave-vector dependence of spin and charge susceptibilities. In 
Sec. we conclude, with remarks on the perspective of the new scheme concerning its applicability in quantum 
many-body physics. 



2 



k 3 ,s 



ko,s' 



kd, s' 



FIG. 1: The coupling function V(ki, k2, fe). Wavevector k$ is fixed by wavevector conservation on the lattice. 



II. THE METHOD 



We start with a grand-canonical Hamiltonian of the form 

H = Vf(fc) c t c r + V V(k,k',k + q)ct ct cv, ,c t 

Z-^ sv ; k,s k,s oat Z—i V ' ' *J k+a,s k'-n,s' k ' ,s' fc,, 



(1) 



The eg ^-operators annihilate electrons with wavevector k and spin projection s z = ±1/2. /i is the chemical potential, 

£(fc) = e(k) — \i the free dispersion of electrons hopping on a 2D square lattice measured with respect to the free 
Fermi surface, and V(k, fc', q) defines the spin-rotationally invariant and wavevector-conserving interaction between 
electrons. For the Hubbard onsite-repulsion we have V(k, k', q) = U > 0. 

The interaction flow scheme is derived using the functional rcnormalization group formalism for one-particle ir- 
reducible (1PI) vertex functionsAi^ It describes the change of 1PI vertex functions exactly upon variation of a 
parameter £ in the quadratic part of the action Qf . Writing Grassmann fields for the fermions and using a combined 
index fc = {iko, k) for Matsubara frequencies fco and wave vectors fc, with the Hamiltonian Q the full quadratic part 
of the action reads 



Cfc 



Ck ■ 



(2) 



In the momentum-shell RG Qi is defined by supplementing Q with a cutoff function depending on the running 
infrared cutoff I. The temperature-flow scheme^ uses the temperature T appearing in front of the frequency sum 
and in the Matsubara frequencies as flow parameter. The change of all vertex functions with I is determined by an 
infinite hierarchy of equations, obtained from the ^-dependent Legendre transform 
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of the generating functional for connected, non-amputated m-point correlation functions W^ r ' 
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(4) 



with <S(c, c, C) being the standard action used to generate correlation functions. 

We condense the notation and use a combined index p — (iui n , fc, s), moreover J dp — T^2 n s j mJj2 ■ The infinite 
hierarchy of equations is truncated in the usual way by setting the six-point function to zero for all t. The equations 
describing the evolution of the self energy Yii(p) and the coupling function Vi(pi,p2,P3) with p± = p\ +P2 — P3 (from 
which the fully antisymmetric four-point vertex can be reconstructed^, see also Fig. read 



— £^0) = / dp'St(p') [Ve(p,p',p')-2Ve(p,p',p)] 
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FIG. 2: The particle-particle and particle-hole diagrams contributing to the one-loop RG equation for the interaction vertex. 
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with 



jgVl(pi,P2,P3) = Tpp.l + Tp H ,l + Tp r H,i 



dpVt(px,P2,p) L(p, -p + Pl +p 2 )Vi(p, -p+Pl +P2,P3) 



TphApi,V2\Pz,P4) = 



dp 



^Vi{px,p,p3) L(p,p + pi -p a ) Vi{p + pi -pz,P2,p) 
+V e (pi,p,p + pi -p 3 )L(p,p + p 1 -pz)Vi{p + pi -P3,P2,P) 
+Ve(p 1 ,p,p 3 ) L(p,p + p x ~p 3 ) V e (p 2 ,p + Pi ~P3,p) 

Tp H £{p\-,P2\P3-,P±) = 

dpVz(pi,p + p2 -p 3 ,p)L(p,p + p 2 -P3)Ve(p,p 2 ,P3) 
In these equations, the product of the two internal lines in the one-loop diagrams is 

L(p,p') = S t (p)WP(p') + wP(p)S e ( P ') 
with the so-called single-scale propagator 
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The one- loop diagrams corresponding to the terms |Q, iJSJ and © arc shown in Fig.|3 

To obtain the desired interaction flow scheme we first multiply Q with a scale factor 1/g and split it in two, yielding 
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(12) 
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g will be the flow parameter. We can absorb the factor 1/g in rescaled fields c, c defined as c = g x l 2 c. As a 
consequence the quartic interaction term 

= — ^ V(k,k',k + q)c k+q ^c k >^ q ^>c k ' :S 'Ck,s (13) 
picks up an extra factor g 2 when written in terms of the new fields: 

V 9 4) = Jn ^ g 2 V(k,k',k + q)h +q ,sh'- q ,si£ k >,si~c k , s . (14) 

k ,k ,k-\-q 

Higher order interactions ]/( 2 ™ 1 ) are rescaled according to Vg 2m ^ — g m V^ 2m \ meaning that changing the scale factor 
1/g in Q g corresponds to changing the strength of the bare interactions. The rescaled fermions c, c describe a system 
with a bare interaction strength g 2 V . In particular, we can start at an infinitely small g, i.e. at infinitely small bare 
interaction, and use the RG formalism to integrate up to the desired bare interaction strength, reached at g — 1. We 
can also stop the flow at any other value of g, with the functions gT, and g 2 V g (k 1 k', k + q) being the self energy and 
interacting vertex function for the bare interaction g 2 V(k,k' ,k + q). Most importantly for practical purposes, the 
truncation of the RG hierarchy and the subsequent neglect of self energy feedback ( see below ) does not change this 
property. 

We call this the interaction flow (IF) scheme. Compared to straight-forward perturbation theory, the advantage 
of the functional RG machinery in the one-loop approximation is that it intrinsically sums up all one-loop diagrams 
contributing in perturbation theory. Therefore, all kinds of ladder and bubble summations, as well as the corresponding 
vertex corrections are included. 

We note that, opposed to the cutoff or T-flow RG schemes, singularities on the right hand side of the flow equation 
are not regularized by the flow parameter. For example, the particle-particle bubble appearing in the perturbation 
expansion will typically diverge logarithmically when T — > 0. Thus, the IF scheme has to be performed at finite 
temperature when the individual one-loop diagrams are bounded. The strength of the interaction is increased 
continuously, and in the course of the interaction flow potential singularities are approached from below. Fig. |3| 
illustrates how the various methods detect perturbative singularities from different directions in parameter space. We 
further remark that the IF scheme does not correspond to viewing a system on different length scales. This implies 
that, unlike in the cutoff RG approach, phase transitions cannot be directly related to fixed points of the flow, where 
the action becomes invariant upon rescaling of fields and spatial coordinates. 

Propagators: The one- loop diagrams consist of a single-scale propagator S g (p) and full two-point Green's functions 

(2) 

Wg (p) . When we neglect self energy corrections they are given by 



WP( P )=g 



-iw„ 



(15) 



and 



S g {p) = -W^(p)Q g W^{p) = 



(16) 



We have arrived at a compcllingly simple prescription for the one-loop flow: we can compute the diagrams just as in 
perturbation theory and only need to multiply them with a scale factor g m ~ l . This shows a striking advantage of the 
IF scheme concerning the computational effort if self energy corrections are neglected: the one-loop corrections only 
depend on the flow parameter through a simple multiplicative factor <?" l_1 , and therefore the corresponding diagrams 
have to be calculated only once in the beginning of the flow. This is in contrast to the cutoff or T-flow schemes, 
where, due to the absence of scale-invariance, the one-loop terms have to be recalculated in every RG step. 

We note another crucial difference between the IF scheme and the cutoff or T-flow RG: in the latter the single-scale 
propagator receives its dominant contributions from a limited shell in the Brillouin zone, either defined by the 
scale-derivative of the cutoff function, or by the T-derivative of the Fermi function. In the course of the flow this 
shell contracts around the Fermi surface, such that at low scales or temperatures only quasiparticles near the Fermi 
surface can contribute. For weakly coupled or Fermi liquid-like systems, these are expected to be sufficiently well 
defined. Hence, neglecting self energy effects such as finite lifetime can be a good approximation^ In the IF scheme, 
however, S g (p) is only weakly confined to a shell around the Fermi surface by the band energy in the denominator. 
This may lead to an overestimate of the influence of Brillouin zone regions away from the Fermi surface and should 
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FIG. 3: Depiction of the different functional approaches in the parameter space spanned by interaction strength g, temperature 
T and infrared cutoff A. The surface represents the critical manifold below which the perturbation theory diverges. 

be corrected in future improvements of the method by inclusion of lifetime effects. 

Simple example: To see how infrared singularities are treated in the IF scheme, consider a simple Cooper instability 
for a local attraction Vo < between electrons. This coupling constant will be subject to one-loop corrections when 
the scale factor g is increased, starting from an infinitely small initial value. We restrict the flow to the particle-particle 
channel for electrons with incoming momenta k and —k. In this case the diagram corresponding to Tpp t i(k, —k, k! , — k') 
leads to a flow of the scattering vertex V g (k, — k, k') = V p given by 



— V p = -V p 
dg 9 9 



dp 



S g ( P )W^(-p) + W^(p)S g (-p) 



-V q p 2g XPP (0) 



(17) 
(18) 



where xpp(O) = / u' 1 +e(f)' i * s ^ ne particle-particle bubble for zero total momentum and frequency. This equation 
can be integrated from to g, yielding 



9 2 V 



9 2 V 



i + g 2 v oX pp{0) 



(19) 



The Cooper instability occurs for g 2 Vo = — 1/xpp(0), corresponding to the usual criterion at g = 1. This illustrates 
how the instability is approached by increasing the bare coupling g 2 Vo. The critical bare coupling depends on 
temperature through the temperature dependence of the one-loop diagrams, which in the Cooper case is xpp(O). 
Since xpp(O) diverges logarithmically for T — > 0, the critical coupling becomes infinitely small for T — > 0. 

All other typical infrared instabilities due to summations of bubble or ladders are contained in the IF scheme in a 
similar way. In particular, particlc-hole instabilities at small wave vectors q, such as ferromagnetism or Pomeranchuk 
Fermi surface deformations-, are treated on equal footing with Cooper and large-g instabilities. Thus, the IF scheme 
combines advantages of the temperature-flow RG with a very lucid formalism closely related to standard perturbation 
theory. 

Effective coupling constants in charge, spin and pairing channels: It is useful to define certain classes of coupling 
functions, which drive certain susceptibilities through one-loop corrections. If any of these classes diverges, it signals 
a potential instability in the corresponding channel. We define effective couplings 



V c (k,k',q) = 2V(k,k',k + q) - V{k',k,k + q) 



(20) 
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in the charge channel, 

V s (k,k',q) = -V{k',k,k + q) (21) 

in the spin channel, and 

V sc {k, k', q) = V{k, -k + q, k') (22) 

in the pairing channel. The flow of these quantities can be analyzed, and their wavevector dependence contains 
information on the dominant tendencies of the system. For example, a d x 2_ y 2-w&ve Cooper instability leads to a 
d x 2_ y 2 dependence of Vsc(k, k' , q — 0) on k and fc', not necessarily of a separable form. 

We have dropped the frequency by writing k instead of k, in anticipation of the following applications in which we 
neglect the frequency dependence of the interaction vertex. 

Coupling to external perturbations, susceptibilities: We introduce fermionic coupling terms to order parameters or 
static external fields corresponding to different symmetry breaking channels. Together with the flow of the interaction 
vertex we compute the one-loop renormalizations of these coupling terms and thereby obtain information on the 
growth of the corresponding fluctuations during the flow. For details about this procedure see e.g. Ref. A 
divergence of one of these external couplings signals a crossover to a region with a strong tendency to ordering. This 
could become an actual finite temperature phase transition if allowed, e.g. when coupling in a third spatial direction 
is added. 

The d-wave pairing field Qd-scik), couples to the rescaled fermions at interaction strength g 2 and wavevector fc via 

with a coupling constant hd- sc (k). The initial condition is taken to be a d x 2_ y 2 form factor hd- sc {k) = 
(cosk x — cos k y ) /V2- As for the other coupling constants defined below, the fc-dependence of hd- sc (k) will evolve 
through corrections of the Cooper coupling Vsc(k,k' ,q = 0) in the course of the flow, and higher harmonics of the 
same representation of the point group as well as a dependence on the distance to the Fermi surface will be generated. 
Notice that we already wrote $ Sy d-sc, related to the source in the IF action via $d-sc = g&d-sc- The susceptibility 
Xd~sc in this channel at interaction strength g 2 is obtained from the susceptibility Xd-sc computed in the IF formalism 
as Xd-sc = Xd-sc/ g 2 - One can check this relation by considering the non-interacting case, where g merely rescales the 
free fermionic Green's function. 

We also define spin-density and charge-density wave fields & s (k,q) and & c (k,q) with point group symmetry 
They are renormalized by one-loop corrections involving the spin and charge couplings V s {k,k',q) and V^fc, fc', (7), 
respectively. Due to spin rotation invariance it is sufficient to consider spin fields in the spin quantization direction, 
coupling via 

¥ ZiS (k, q) h s s (k, q) (% +f)S Cg )S - % +f ,_ s % _ s ) , (24) 

where in the s-wave channel the initial condition is taken to be h s s (k,q) = 1. 

In the 2D Hubbard model the momentum transfer q = Q = (n, tt) between created and annihilated particles plays an 
important role. It corresponds to an alternating field acting on electron spins on n.n. sites, and if the corresponding 
coupling h s s (k,Q) becomes large, this signals a strong tendency towards antiferromagnetic (AF) spin-density wave 
(SDW) order. We also consider couplings in the d-wave channels, in order to look for e.g. d-density wave fluctuations, 
coupling to particle- hole pairs in the charge channel with momentum transfer Q: 

$t(kAY, hd c& E U$,s%s< (25) 

with initial condition chosen as hf(k) — (cosk x — cosk y ) In addition to the (tt, n) channels, the q = 0-channels 

will also be analyzed. The isotropic s-wave part contains information about the compressibility of the Fermi sur- 
face, while higher s-wave, d- and p-wave parts are important when estimating tendencies towards Fermi surface 
deformation^ 
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III. INTERACTING FERMIONS IN ONE DIMENSION 



Here we briefly describe how the IF scheme reproduces some well-known perturbative results for electrons in 
one spatial dimension. Consider the g-ology model with spin-independent interactions where left- and right-moving 
electrons around two Fermi points ±kp are scattered by coupling constants g\ = V(kp, s; — kp, s' — > — kp, s; kp, s'), 
cj2 = V(kp, s; —kp, s' — > kp, s; — kp, s') and 54 = V(kp, s; kp, s' — > kp , s; kp , s'). For densities sufficiently away from 
half filling Umklapp processes g$ — V(kp, s; kp, s' — > —kp, s; —kp, s') can be neglected. One-loop corrections in the 
particle- hole channel with wavevector transfer 2kp and in the zero-total- momentum particle-particle channel diverge 
at low temperatures like \og(W/T), where W denotes the band width. Focusing on these contributions, and writing 
s instead of the scale parameter g in order to avoid a <?-mess, the one-loop IF equations read (using No = l/(irvp) 
for the density of states) 

= -ATo^logf (26) 

— = -- 3l , s2s log- (27) 

<?4 is not renormalized by logarithmic corrections. When we now integrate form zero to the desired bare interaction 
strength reached at s = 1, eq. 1261 immediately yields (with the bare coupling g\ = g\. s =o) 

5l ' S=1 = 1 + Clogf ■ (28) 

This means that at low T <^W a, repulsive g\ is renormalized away by one-loop contributions. The low-temperature 
physics is thus determined by the remaining coupling constant 32 which approaches a constant g% for T — * 0. The 
coupling to external spin- and charge density perturbations with wavevector 2kp, hsuw and /icdw, are given by the 
one-loop corrections 

W 

isdw 2s log — (29) 



d/lSDW 


_ N 


ds 


2 


dhcnw 


No 


ds 


2 



T 

, W 

■s log y 



(30) 



Here we inserted the low temperature fixed-point values for the coupling constants g\ tS and g 2 ,s- This is justified for 
low T, since then g\ and g\ are approached rapidly for small s. For repulsive 171 > and 32 > 0, the SDW coupling 
grows most strongly when the bare interaction strength s is increased. Integrating Eq. 1291 from s = to s = 1, we 
find the T — > power laws for the SDW-coupling and the corresponding susceptibility: 

^sdw(T) cx T- N ^' 2 , xsdw(T) cx T- N °s'2 . (31) 

These powerlaws are in agreement with results obtained via bosonizationi 7 " by expanding the Luttinger liquid param- 
eter K p up to first order in the fixed point coupling g^. Likewise, the leading susceptibilities in the g\-g2 plane are 
reproduced correctly by the IF scheme. 

In repulsive Luttinger liquids the presence of a single impurity leads to universal power-laws and scaling functions in 
the low energy limit. Using a Matsubara frequency cutoff as flow parameter for a spinless fermion model with nearest 
neighbor interaction and a single static impurity, the computation of spectral properties of single-particle excitations 
and oscillations in the density profile induced by the impurity or a boundary was shown to provide remarkably accurate 
results even for moderate coupling strengths in comparison with exact asymptotic and numerical results^ It turns 
out that within the IF approach these effects of a single static impurity in a Luttinger liquid are well captured and 
the respective predictions are confirmed quantitatively. A more detailed comparison is under way. 



IV. APPLICATION TO THE TWO-DIMENSIONAL HUBBARD MODEL 



Here we apply the IF scheme to the two-dimensional t-t' Hubbard model on a square lattice and compare the results 
with those obtained in other RG schemes. The latter have been widely used to identify the leading low-temperature 
instabilities of the Hubbard model. Practically all numerical 2D implemcntations&liiiSiiLifliiii use a patch discretization 
introduced in this context by Zanchi and Schulz^ The coupling function is calculated for two incoming and one 
outgoing wave vector near or on the Fermi surface, and is taken to be constant within elongated patches that extend 
roughly perpendicular to the Fermi surface (see Fig.^J. This approximation is believed to be good, since the leading 
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FIG. 4: Left: iV-patch discretization for the cutoff and T-fiow RG schemes. Right plot: Box discretization used for IF. 
open circles denote the iV 2 points for which the coupling function is calculated. The shading indicates the band energy. 



The 



flow is typically given by the coupling functions connecting particles near the Fermi surface, and because during the 
flow at least one of the internal lines in the one-loop diagrams is restricted to an increasingly narrow neighborhood of 
the Fermi surface. As mentioned in the previous section, the second "contraction" feature is absent in the IF, since 
no cutoff function is used, and no derivatives of the Fermi function occur in the one-loop diagrams. Therefore the 
full BZ contributes to the flow at all values of the flow parameter g, and if the couplings near the FS become large it 
may not be a good approximation to also use them for BZ regions far away from the FS. We are thus forced to use a 
different discretization of the BZ. This actually reduces the degree of approximation caused by the usual projection 
on the Fermi surface. 

We divide the BZ into N 2 squares and calculate the flow for wave vectors in the centers of the boxes, keeping the 
couplings constant within each box. With this prescription we can numerically treat up to iV 2 = 18 2 = 324 boxes. 
Note that this discretization does not correspond to a finite-size system. The internal states in the one-loop diagrams 
are still taken as those of an infinite lattice and can lie arbitrarily close to the Fermi level. 

We focus on band fillings close to van Hove filling and vary the nearest neighbor hopping t' . In the temperature- flow 
RG the picture is the following 10 : For small t' and near half filling the Fermi surface is well nested, and the flow 
to strong coupling at low temperatures is mainly of spin-density wave character. For larger \t'\ and the FS fixed at 
the saddle points, d-wave pairing tendencies become stronger, while SDW tendencies are suppressed. For t' < —0.2t 
d-w&ve pairing dominates over SDW. For t' — ► —t/3 the critical temperature scale for the flow to strong coupling 
drops by several orders of magnitude. For t' < —t/3 it rises again, and the flow to strong coupling is dominated by 
ferromagnetic tendencies. The temperature-flow RG thus suggests a quantum critical point at t' « — 1/3 between a 
d-wave paired and a ferromagnetic ground state. 

At low temperatures, the IF results show (roughly) the same regimes when \t'\ is increased, and the FS is fixed at 
the van Hove points. The results for t' — are shown in Figs. |S] and H3 When the interaction strength is increased 
for fixed temperature, we find a critical value g 2 (T) = U C (T) at which the maximal coupling function diverges. This 
value corresponds to the critical interaction strength needed to destroy the weakly coupled state. For the perfectly 
nested case, i.e. t' = and /i = 0, the SDW coupling grows strongest when the bare interaction approaches U C (T). 

In Fig. [7] we compare the characteristic temperature for the spin-density wave instability as a function of the bare 
interaction U for several RG-like schemes with RPA values. In all cases T C (U) is a monotonically increasing function of 
U, due to its mean-field character. The T-flow RG scheme yields the smallest values for T C (U), RPA the highest. The 
momentum-shell 1PI scheme lies slightly below the RPA values, and the IF data are roughly in the middle between 
RPA and T-flow. 

When we move away from perfect nesting by introducing a next-nearest neighbor hopping amplitude t' < 0, 
and change the band filling at low temperature T = O.OOli, there is an increasing tendency towards d x 2_ y 2-waye 
superconducting pairing. For t' — — 0.3t this is shown in Figs. |H1 and In Fig. [HI the band filling is slightly above van 
Hove filling, and different discretizations show consistently that the effective coupling in the d-wave pairing channel 
grows more strongly than the leading magnetic component at wave vector Q — (tt,tt). Fig. [5] shows analogous data 
at the saddle point filling. With increasing patch number the leading divergence in the SDW channel seems to win 
over the pairing instability. 

As already mentioned, ignoring lifetime effects possibly makes us overestimate the influence of excitations away 
from the Fermi surface. In Fig. 1101 we see that the coupling function in the particle-hole channel at wavevector 
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FIG. 5: Left: Flow of the most attractive and most repulsive coupling constant as a function of the bare interaction strength 
g 2 = U for the perfectly nested case t' = 0, fi = 0. Right: Growth of the coupling to static spin perturbations at q = (7r,7r) 
(spin-density wave channel). The two other degenerate lines that rise much slower denote the effective couplings in the d-density 
wave (dDW) and d-wave pairing (dSC) channel. 




FIG. 6: Data for half filling and t' — 0, T — O.OOlt at the critical bare interaction strength U = 0.57t where the couplings 
diverge. Left: (f-dependence of the interacting spin susceptibility. The peak is at q = (tv,tv). Right: ^-dependence of the 
interacting charge susceptibility. Note the difference in the scale. 



transfer q w (ir,ir), i.e. the SDW channel, is strongly enhanced in all BZ regions, and is largest at the top and 
bottom of the band. It is likely that in these regions with large band energy the lifetimes of quasiparticles will 
be short, and the flow without self energy corrections may overestimate the SDW tendencies. In contrast to this, 
the pairing interactions become large only in the vicinity of the Fermi surface. Only a part of this fc-dependence 
is due to the d-wave form factor. Due to the discretization scheme we might actually underestimate the growth 
of the pairing correlations, since the coupling function cannot be calculated arbitrarily close to the Fermi surface 
using only 18x18 boxes. Hence, we cannot draw definite conclusions about the leading divergence directly at van 
Hove filling. Slightly away from van Hove filling d-wave tendencies grow more rapidly than SDW tendencies. This 
trend is consistently reproduced for all discretizations we used. Keeping in mind that we presumably overestimate 
the SDW channel and underestimate the pairing tendencies, we find that the dominance of d-wave pairing over 
SDW (as shown in Figs. 181 and llUfl is likely to be more robust than suggested by the flow without self energy corrections. 

In Fig. ^2 we show data obtained within the cutoff RG for the one-particle irreducible vertex functions from Ref. 
Ilii using the same discretization scheme as in the IF algorithm. There is good agreement with the IF data in Fig. 
ITUl We notice, however, that the SDW channel is a little weaker, and the effective couplings to external SDW fields 
are now peaked (albeit very mildly) around the FS. Again, this difference occurs because in the cutoff scheme the 
modes contributing to the flow are restricted to narrow shells around the Fermi surface, meaning that at small scales 
only those couplings involving modes near the Fermi surface grow most strongly. 

Very close to van Hove filling the overall behavior of the IF flow resembles the so-called saddle point regime of Refs. 
[H andB The dominant scattering processes near the Fermi surface are those connecting the saddle points, and 
SDW and d-wave pairing susceptibilities both become large. Furthermore, the charge compressibility is suppressed, 



10 




U/t (t/U) 1/2 



FIG. 7: Comparison of the characteristic temperatures for the SDW instability at t' = 0, n = for several perturbative RG-like 
schemes. The right plot shows the same data with different a;-axis. Squares: RPA (360 x 360 lattice in the BZ), circles: T-flow 
(48 patches around the FS), triangles: momentum-shell- RG (lPI-scheme, 48 patches around the FS), crosses: IF for 144 patches 
in the BZ. 
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FIG. 8: Results for the d-wave Cooper instability at fj. = — Lit, t' = — 0.3t and T = O.OOlt. Left: Flow of effective coupling in 
ci-wave pairing (solid line) and AF SDW channel (dashed line) for a 18x18 discretization. Right: Ratio of these two effective 
coupling constants for different discretizations: solid line 18x18, dashed-dotted line 16x16, dashed line 14x14, dotted line 
12x12. 

predominantly in the Brillouin zone regions around (jr,0) and (0, 7r), similar as in Ref.E This can be seen from 
the flow of the /c-dependent couplings to external static charge perturbations at long wave lengths, h c [k, q— ► 0) (see 
Fig. I12|) . A similar picture is found at half filling. Although at weak coupling these tendencies are rather small, they 
connect qualitatively well with much stronger effects found by Otsuka et in Quantum Monte Carlo simulations 
for U = At and T = 0.2t. 

Upon increasing the temperature in this d-wave regime, the instability is pushed to larger bare interactions. For 
t' = — 0.25i and (n) = 0.9 the <i-wave instability occurs for temperatures below T = 0.0035i at U < 2.75t. This is 
in reasonable agreement with temperature-flow RG results (with a completely different discretization of the Brillouin 
zone), where the flow diverges at this temperature for U ~ 2.6t. For higher temperatures, the leading divergence 
appears to occur in the SDW channel at wave vector (tt, tt) again, and the interaction needed for the instability barely 
increases with rising T. This indicates that spin fluctuations live on a larger energy scale than pairing tendencies. 
A similar change in the character of the flow to strong coupling with increasing temperature has been observed in 
cutoff RG schemes^ A conservative interpretation suggests the presence of strong spin fluctuations in the normal 
state above the superconducting ground state (when 3D coupling is present such that long range order is possible). 
An improved treatment including lifetime effects will give further insights concerning the strength of these effects. A 
general trend within the present approximations is that in the IF scheme the SDW tendencies seem stronger than in 
other RG approaches. 

When \t'\ is increased further, the critical interaction strength for the flow to strong coupling increases, and for 
fixed U ~ 3t the couplings remain bounded down to very low temperatures. For t' < — 0.36t we find that the main 
instability is in the ferromagnetic channel. The change in the dominant spin correlations can also be seen in Fig. 
1141 where we plot the spin susceptibility for different t' at van Hove filling along a trajectory in (f-space. The results 
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FIG. 9: The same for van Hove filling at fi = —1.2t, t' = — 0.3t and T = O.OOlt. The (ft, tt) spin channel may win for large 
patch numbers. 




FIG. 10: IF data for T = O.OOlt and (i = -0.8t, t' = -0.25* «n) = 0.9). Upper left plot: Flow of the average coupling to 
external d I 2_ !/ 2-pairing fields, hd-sc (solid line), and SDW fields, hl(q = (rr, 7r)) (dashed line), corresponding to the dominant 

tendencies in this case. Upper right: d-wave coupling hd- BC (k,q — 0) (solid lines with crosses) for k = (k x ,0) and k — (k x ,ir) 
at the critical interaction strength. SDW coupling h s s (k,q — (ir,ir)) (dashed lines) for k = (k x ,0) and k = (k x ,n). Lower 
left/right: d-wave and SDW couplings for k varying through the full Brillouin zone. The scalebar depicts the strength of the 
coupling. 

agree qualitatively well with those of the temperature-flow scheme, which finds the onset of the ferromagnetic regime 
at t' « -0.33* for U = 3t. 



V. CONCLUSIONS AND OUTLOOK 

We have presented a perturbative method to study interacting fermionic systems, the so-called interaction flow 
(IF) scheme. Based on the functional renormalization group formalism^ it sums up all one-loop corrections in an 
unbiased manner using the bare interaction strength as a continuous flow parameter. We have applied the IF scheme 
to one-dimensional fermions and the two-dimensional Hubbard model and compared the results with those obtained 
by cutoff RG and temperature-flow schemes. There is ample qualitative and in some cases reasonable quantitative 
agreement, even when very different Brillouin zone discretizations are compared. 
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FIG. 11: Cutoff RG data, again for T = O.OOlt and /j, = -0.8t, t' = -0.25*, [7 = 2.5*. Upper left plot: Flow vs. infrared 
cutoff A of the average coupling to external d x i_ y i -pairing fields, /id_ sc (solid line), and SDW fields, h s s (q = (n,iv)) (dashed 
line), corresponding to the dominant tendencies in this case. Upper right: d-wave coupling hd- S c{k,q = 0) (solid lines with 
crosses) for k = (k x ,0) and k = (k x ,n) at the critical interaction strength . SDW coupling h a s (k,q = (dashed lines) for 

k — (k x ,0) and k = (k x ,n). Lower left/right: d-wave and SDW couplings for k varying through the full Brillouin zone. The 
scalebar depicts the strength of the coupling. 




\ /a U / 1 

FIG. 12: IF Data for T = O.OOlt and van Hove filling \i = -t, t' = -0.25t. The left plot shows the coupling to external 
long- wavelength charge perturbations h c (k, q — » 0) for k varying through a quarter of the Brillouin zone. The tick lines in the 
right plot shows the flow of h c (k, q — * 0) at two points k. The solid line is for a point close to a saddle point k « (ti, 0) and the 
dashed line for k ~ (1.3, 1.3) in the BZ diagonal. The thin lines show analogous data for q = (n,n). 

In agreement with other RG approaches the IF scheme yields an extended region in parameter space of the 2D 
i-<'-Hubbard model where d-wave superconductivity is the dominant low-temperature instability, though it is less 
prominent, and the tendencies towards SDW ordering appear to be somewhat stronger. This holds in particular 
for van Hove filling and t 1 < —0.2t, where, opposed to temperature and cutoff flow approaches, the IF scheme does 
not find a d-wave dominated instability. We believe that this difference is due to the fact that in the IF scheme all 
parts of the BZ contribute to the flow to strong coupling. Neglecting self energy corrections may therefore lead to an 
overestimate of the SDW channel, whose couplings grow over the whole Brillouin zone. In contrast to this, interactions 
driving d-wave superconductivity become large only in the vicinity of the Fermi surface. This problem does not occur 
in cutoff RG or temperature-flow schemes. In these approaches only modes in a contracting shell around the Fermi 
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FIG. 13: Summary of the flow at van Hove filling as a function of i! for T = O.OOlt and a 12 x 12 discretization. U c denotes 
the critical value of bare interaction strength U at which some components of the coupling function exceed the perturbative 
range when U is increased from below. With increasing \t'\, for —t' < 0.2t the instability is of AF-SDW type, then it becomes 
a d-wave like Cooper instability. Note that for larger system sizes the boundaries of the SDW-d-wave crossover shifts (see 
discussion in the text). For — t 1 > 0.36i the dominant instability is ferromagnetic. 
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FIG. 14: Interacting spin susceptibility Xs(o) at van Hove filling and T = 0.0014 for different t' and for the critical interaction 
strength where the couplings diverge, using a 12 x 12 discretization. For t' > — 0.35t, Xs{o) is peaked at q = (tt, tt) = M and 
for t' < -0.35t it is peaked at q = (0, 0) = O. 

surface contribute in the final stage of the flow. This way, the IF scheme sheds some light on the inner workings of 
the previous RG approaches. The inclusion of lifetime effects in the IF scheme will clarify these issues further. 

Furthermore, for large absolute values of the nearest neighbor hopping t' the leading instability is in the ferromag- 
netic channel. This confirms results obtained within the temperature-flow RG and other perturbative approaches^ 

As new results we have presented detailed wave vector dependences of interactions and susceptibilities over the 
whole Brillouin zone. 

We emphasize that we consider the interaction flow approach an additional tool for the analysis of interacting 
many-fermion systems. Its concept is as simple as perturbation theory, yet it captures many of the non-trivial 
effects crucial for the understanding of interacting electron systems. Most importantly the IF scheme includes the 
coupling between the various one- loop channels. This provides an alternative practicable method to detect and gauge 
infrared instabilities of fermions in low dimensions. Having several viable techniques at hand may prove useful, as the 
character of different approximations can be analyzed more clearly. " True" renormalization group methods like cutoff 
RG schemes appear to have a broader range of applicability, and seem to be the wider concept. For example, RG 
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methods establish a relationship between systems on different length or energy scales in a quite general way, while 
the IF flow is by construction restricted to a perturbative regime. 

Regarding future applications and extensions of the IF scheme we note the possibility of comparing the results 
with Quantum Monte Carlo (QMC) data in the 2D Hubbard model and on ladder systems. BuhrUi gives detailed 
results for the irreducible vertex function at U — At in the two-leg ladder and the 2D Hubbard model, including wave 
vector dependence. With the current IF scheme we can easily achieve a comparable fc-space resolution. Hence, it 
will be interesting to see how the perturbative method compares with QMC. This can give new insights on how the 
QMC data extrapolate to lower temperatures, that cannot be reached by QMC due to the sign problem. Moreover, 
we believe that the inclusion of Fermi surface renormalizations and distortions, as well as lifetime effects are better 
tractable in the IF scheme, since complications due to a flowing cutoff do not occur. The example of an impurity in a 
Luttinger liquid described briefly at the end of Sec. IIIII shows that the inclusion of self energy effects in the IF scheme 
can lead to excellent results. 

All in all, we hope that the IF scheme will open new possibilities in the understanding of interacting many-body 
systems. 
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